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Abstract 

Explicit expressions for inverse of Young's modulus E(n), inverse of shear modulus G(n, m), and 
Poissons ratio z/(n, m) for cubic media are considered. All these characteristics of elastic media 
depend on components Sn, Su and S44 of the compliance tensor S, and on direction cosines 
of mutually perpendicular vectors n and m with fourfold symmetry axes. These characteristics 
are studied for all mechanically stable cubic materials for vectors n belonging to the irreducible 
body angle subtended by three cubic high symmetry directions [001], [111], and [110]. Regions 
in the stability triangle of in which cubic elastic materials are completely auxetic, nonauxetic, 
and auxetic are established. Several intermediate-valence compounds belonging to the region of 
complete auxecity are indicated. The extreme properties of E~ x , G^ 1 and v established by Hayes 
and Shuvalov are confirmed. 
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I. INTRODUCTION 



Mechanical properties of elastic media are described by the bulk, Young's and the shear 
modules. Poisson's ratio is another important mechanical charactersitics of such materials. 
With the exception of bulk modulus, generally, these characteristics are anisotropic, i.e. they 
depend on direction of load and direction in which the response of a material is measured. 
The Poisson ratio is especially interesting, because it can be negative. Materials which 
exhibit a negative Poisson's ratio are referred to as auxetic media [1, 2]. Auxetics respond 
to imposed uniaxial tension with lateral extension in place of expected contraction. Lakes 
described the synthesis of an actual auxetic material and proposed a simple mechanism 
underlying the negative Poisson's ratio [3]. Such media can find interesting applications in 
future technologies. Ting and Barnett introduced classification of the auxetic behavior of 
anisotropic linear media [4]. 

Among the linear elastic media, a special attention has been devoted to cubic media 
[1, 2]. Turley and Sines used the technique of rotations by Euler angles [5], and considered 
anisotropy of Young's modulus, shear modulus, and Poisson's ratio in cubic materials. Using 
the technique developed by Turley and Sines, Gunton and Saunders studied stability limits 
on the Poisson's ratio, and considered martensitic transformations [6]. Milstein [7, 8], Huang 
et al. [9] analysed responses of ideal cubic crystals to uniaxial loadings. 

Poisson's ratio for cubic materials was studied by Milstein and Huang [10], Jain and 
Verma [11], Guo and Goddard [12], Erdacos and Ren [13], and Ting and Barnett [4]. With 
the exception of the last paper, the authors were looking for directions for which Poisson's 
ratio is negative. In the opinion of Baughman et al, negative Poissons ratios are a common 
feature of cubic media [1]. They refer specialy to cubic metal single crystals. Our results 
support this opinion [14]. 

An attempt to study more global properties of cubic auxetics was made by Tokmakova 
[15]. To describe elastic properties of cubic media, she used the familiar elastic parameters 
C11/C44 and C12/C44. The purpose of our paper is to study anisotropy of all three mechanical 
characteristics of all mechanically stable elastic media. To achieve this goal, we shall use 
different parameters, which were introduced by Every [16] and which were used by us for 
the unified description of elastic, acoustic and transport properties of cubic media [17-20]. 
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II. ELASTIC PROPERTIES OF CUBIC MEDIA 



Consider homogeneous mechanically stable anisotropic elastic material. It is characterized 
by the components S^i {ij, k, I — 1, 2, 3) of the compliance tensor S, or by the components 
Cijki (ij, k,l = 1, 2, 3) of the stiffness tensor C. Both these tensors have the structure 
[[V 2 ] 2 ] [21], i.e. they have the complete Voigt's symmetry, e.g. C ijk i = C jiM = C Hyij . The 
compliance and stiffness tensors are mutually inverse, i.e. 

S = C" 1 . (1) 

After Walpole [22], we denote degenerate eigenvalues of C by cj, cl, and cm- One has 
[23] 

cj = Cn + 2Ci2, cl = 2C44, cm = Cn — Ci 2 , 

s.j = c~j\ s L = cl 1 , s M = cm 1 . (2) 

Similarly sj, sl, and sm are eigenvalues of S. 

The eigenvalues sj, sl, and sm can be written in terms of S u , Su and 544, namely, 

sj = S n + 2Si2, s L = 544/2, s M = S n - S 12 . (3) 

Using these relations and the second line of Eqs. (2), we express the matrix elements Suv 
bjC uv (U,V = 1,2,4) 

Sll = , 5*12 = , 5*44 = — — . (4) 

CjCm CjCm C/44 

These expressions are in agreement with the familiar results [24]. 

One may use eigenvalues cj, cl, and cm in place of three elastic constants Cn, C12, and 

C44, but the more reasonable choice is the use of three parameters introduced by Every [16], 

namely, Si = (Cn + 2C 44 ) and two dimensionless parameters 

(Cn — C44) (Cn — C12 — 2C 44 ) 

S2 = ,s 3 = . (5) 

Si Si 

The inverse relations hold [20] 

Cn = si (2s 2 + 1) /3, C 12 = si (4s 2 - 3s 3 - 1) /3, C 44 = s x (1 - s 2 ) /3. (6) 

The elements of C and eigenvalues cu (U = J,L,M) are proportional to s\. From Eq. (1), 
it follows that Sij oc s^ 1 . Therefore, we introduce the dimensionless components of the 
compliance tensor S'^ = SiSij. 



The mechanical stability is guaranteed when all introduced eigenvalues are positive. This 
condition yields familiar inequalities [25] 

(A u + 2A 12 ) > 0, (A n - A 12 ) > 0, A u > 0, (A = C, S). (7) 

In terms of Every's parameters, cubic material is stable, if s± > and s 2 , S3 belong to the 
stability triangle (ST for short) in the plane (s 2 , S3) (cf. Fig. 1) 

(1 + 2s 2 ) > |4s 2 - 3s 3 - 1| , (10s 2 - 6s 3 ) > 1, s 2 < 1. (8) 



Trojstab.png 



FIG. 1: Triangle abc of stability (ST) of cubic media. Points (s\, s 2 , S3) characterizing mechanically 
stable cubic materials belong to the prism (53 > 0), the base of which is the stability triangle. 

Introduce two parameters of elastic anisotropy 

Xc = (c M - c L ) = (Cn - C12 - 2C 44 ) = S1S3, 
Xs = si (s M - s L ) = -s 1 [ = = . (9) 



.c M c l j c l c m c l c m 
The parameter Xs is dimensionless. 

For s 3 = both Xc, and Xs vanish, hence the anisotropic parts od E -1 , v and G^ 1 vanish 
(cf. Sect. IV). In this case the obtained expressions for mechanical characteristics are the 
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same as for isotropic materials (cf.[23]). Therefore the interval dc (1/10 < s 2 < 1, S3 = 0), 
indicated by the line lying below the s 2 -axis, is the interval of isotropic materials. We shall 
demonstrate later that the isotropic auxetic materials of cubic symmetry are located on the 
interval 1/10 < s 2 < 1/4 (s 3 = 0) indicated by the interval lying above the s 2 -axis. 



III. MECHANICAL CHARACTERISTICS OF ELASTIC MEDIA 

The response of elastic materials to hydrostatic pressure is characterized by the bulk 
modulus K [24]. For cubic elastic media 



K=[3 (S n + 2S 12 )]- 1 = (3SJ)- 1 = (C u + 2C 12 ) /3. 



(10) 



In Fig. 2 we present the dependence of dimensionless bulk modulus Kj s\ on the parameters 

S 2 , s 3 



bulk-modulus . png 



FIG. 2: Contour map of dimensionless bulk modulus K (s 2 , s 3 ) /s\ 



We shall characterize the response of the elastic material to uniaxial tension in a direction 
n (nn = 1) by the dimensionless Young's modulus e(n) = E(ja)/s\ 



e 1 (n) = (n <g> n) • S' ■ (n ® n) = ninjS' ijkl n k ni. 



(11) 



Repetition of a suffix in a product of tensors or in a single tensor implies the usual summation 
with the respect to that suffix over the values 1,2,3. 

Similarly, the response to the shear stress will be characterized by the dimensionless shear 
modulus g(n, m) = G(n, m)/si, (mm = 1, nm = 0) 

1 



4g(m, n) 



(m <g> n) • S' ■ (m <g> n) = m, i n j S' ijkl m k ni. 



(12) 



Poisson's ratio z/(n, m) is the ratio of the lateral contraction along direction m to the lon- 
gitudinal extension along direction n perpendicular to m 

< \ < \< ^ \ q> ( ^ \ rnimjS' ijkl n k ni 
v(m, n) = — e(n)(m <g) m) • S -(11(8)11) = — — . (13) 

Using the basis of tensors [[V 2 ] 2 ] constructed for cubic media by Walpole [22], we obtained 
the explicit expressions for E(n), [4G(m, n)] -1 , and [— i/(m, n)/E(n)} for media of high and 
middle symmetry classes [23]. For cubic materials, we obtained the familiar Nye [24] and 
Ting and Barnett formulae [4] 

E-^n) = [(sj - s M ) /3 + s L ] + (s M - s L )p(n) = 
S u - 2 (S n - S 12 - ^44/2) T(n) = S u [1 - 3T(n)] + (S u + 2S 12 + S, 4 ) T(n), (14) 

where p (n) = Y^i= a n \i an( i ^( n ) — ( n l n b + n l n l + n c n a)- We also denoted the fourfold 
axes of cubic media by a, b, and c and, for example, m a = (ma). 
Similarly, for z/(m, n) and G(m, n) we get 

- ^^tV = & ~ s m) /3 + % - s L ) P m, n), (15) 
[4G(m, n)]" 1 = s L /2 + {s M - s L ) P(m, n). (16) 

The function P (m, n) was introduced by Every [16] (cf. also Ting and Barnett [4]) 

c 

P (m, n) = P (n, m) = ^ (m^) 2 . (17) 

i=a 

Introduced dimensionless mechanical characteristics of elastic materials are functions of five 
variables. This fact complicates studies of their properties. 

The functions p(n) and P(m, n) obey the familar inequalities 

1/3 =p((lll» < p(n) < 1 =p((001)) , (18) 
= P ((100) , (001)) < P (m, n) < 1/2 = P ((-110) , (110)) . (19) 

Note that the first term of expression defining P~ 1 (n) and [—u(m,n)/E(n)] is isotropic, 
whereas the second one is anisotropic. Shear modulus has the same structure. 

In the Appendix we write down the explicit formulas for three considered anisotropic 
mechanical characteristics in terms of components of both stiffness and compliance tensors 
and consider the limit of isotropic elastic materials. 
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IV. ANISOTROPY OF MECHANICAL CHARACTERISTICS OF ELASTIC MA- 
TERIALS 



Using Eq. (3) expressions for l/e(n), l/g(m, n) and z/(m, n) can be cast to the form 
used by Hayes and Shuvalov [26] 

e-\n) = S[ 2 + S'j2 + X sj:nt, (20) 

i=i 

g-\ m , n) = S' u + A Xs £ m]n\, (21) 

i=l 

t x T (m, n; s 2 , S3) , ooA 

1/(m ' n)= " ^ + gw2 +x ,E^ n r (22) 

where T (m, n; s 2 , s 3 ) = (s 2 , S3) + (s 2 , s 3 ) £Li m 2 n 2 . 

The dimensionless quantities - Young's modulus, shear modulus and Poisson's ratio - are 
functions of angles and Every's parameters s 2 , S3. Studying their angular dependence, we 
may confine ourselves to the body angle Q c subtended by three vectors of high symmetry 
ni = [001], n 2 = [111], n 3 = [110]. An arbitrary unit vector n belonging to f2 c is a linear 
combination of these three vectors 

n = ani + /3n 2 + 7n 3 (nn = 1). (23) 

All unit vectors belonging to the body angle 4-7T can be obtained as a result of the application 
of cubic symmetry operations to vectors belonging to Q c . 

Contour maps of e (n, s 2 , s 3 ), g (n, s 2 , s 3 ), and v (m, n, s 2 , s 3 ) for vectors r^, n 2 , and n 3 
are shown in Fig. 3. 



V. INEQUALITIES FOR ELASTIC CHARACTERISTICS OF CUBIC ELASTIC 
MEDIA 

Properties of mechanical characteristics of elastic media depend on the sign of Xs- I n the 
lower part of ST (where s 3 < 0), Xs > 0, whereas on the upper part of ST, (s 3 > 0) Xs < 0. 

Hayes and Shuvalov studied extreme properties of mechanical characteristics of cubic 
elastic media. Let us summarize their results [26]. If Xs > 0, then at each point of ST the 
function e (n, s 2 , s 3 ) reaches a local maximum e^ x for n = (111), while for n =< 001 > at 
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e-nl-n3 .png 



FIG. 3: Anisotropy of mechanical characteristics of elastic media: contour maps of e (11,52,53), 
g (n, s 2 , S3), and v (n, s 2 , S3) for n = ni = (001), n = n 2 = (111) and n = 113 = (110) 



each point of ST it attains a local minimum e^} n 



,(+) 



(s 2 ,s 3 ) 



- = if> (82, s 3 ) (n=(lll», 



+ 2S' l2 + 5*44 



,(+) 



V«2, S3) 



1 

°11 



^ = F 2 (e) ( S2 , S3 ) (n=(001)) 



If Xs > 0, at each point of ST, value of e for n £ f2 c obeys the inequality e^ n < eW(n) < 
e m2x- Reliefs defined by values of functions (s 2 , S3) and F 2 e ^ (s 2 , S3) intersect along the 
isotropy interval (S3 = 0; 1/10 < s 2 < 1), where neither e nor g (nor 1/) do not depend on 
angles (cf. Fig. 4). 



If Xs < 0, then 



41 (s 2 , s 3 ) = F 2 w (s 2 , s 3 ), (n = (001)) , 



.(e) 



eL ) n(s 2 ,s 3 )=F 1 (e) ( S2 , S3 ),(n = (lll)). 

Relief for local minima of e (n, s 2 ,s 3 ) is obtained with the help of function F^J n = 
( s 2, S3) 0(x s ) + -^i^ (s 2 , S3) 0(— Xs), where 9(x) is the Heaviside step function. 
Consider extreme properties of dimensionless shear modulus. If Xs > 0, then 



9; 



(+) 



= ir = F i ] (*2, s 3 ) (m = (100) , n = (011) ; n = (100) , m = (011)) 



j 44 



Fe-on-isotropy-line . png 



FIG. 4: Contour map of the surface of local maxima of dimensionless Young's modulus defined by 
the function F%L = F[ e) (s 2 , s 3 ) 9( Xs ) + F 2 (e) (s 2 , s 3 ) 6(- Xs ). 

£1 = 0(qf 1 qf ; = F^ 9) (s 2 , s 3 ) (m = (110) , n = (-110) ; n = (110) , m = (-110)) . 
If Xs < 0, then 

aLL = F^ 9) (s 2 , s 3 ) , (m = (-110) , n = (110) ; n = (-110) , m = (110)) 
9tl = f1 9) (s 2 , s 3 ) , (m = (100) , n = (011) ; n = (100) , m = (011)) . 

Properties of surface of local maxima of the dimensionless shear modulus are summarized 
in Fig. 5 



Fg-on-isotropy-line . png 



FIG. 5: Contour map of the surface of local maxima of dimensionless shear modulus defined by 
function Fi g) (s 2 , s 3 ) 9{ Xs ) + F 2 [a) (s 2 , s 3 ) 9(- Xs ). 

In the case of v (m, n; s 2 , s 3 ), Hayes and Shuvalov [26] have proved that if X s > 0, at each 
point of the lower part of ST, the Poisson ratio attains a local maximum for m = (011) , n = 
(100), whereas for m = (011) , n = (100) at each point of the upper part of ST it attains a 
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local minimum 



Vma* (s 2 , s 3 ) = -f^ = F„ (s 2 , s 3 ) (m = (Oil) , n = (100) , X s > 0) , 
fmin («2, s 3 ) = ^ (s 2 , s 3 ) , (m = (Oil) , n = (100)) , < 0). 

Contour map of relief of the function F u (s 2 , S3) is depicted in Fig. 6. 



(24) 



F-nu . png 



FIG. 6: Contour map of relief of F v (s 2 , S3). In the upper part of ST, it is related to geometrical 
loci of local maxima of v. In the lower part of ST, the contour map is related to geometrical loci 
of local minima of v. 

In virtue of inequalities (7), the extreme values of z/(m, n) satisfies the inequalities [26] 
(cf. Fig. 6) 



-1 < v min < (xs > 0, S 12 > 0) , 
< v max < 1/2 ( X s < 0, Sis < 0) . 



VI. PROPERTIES OF POISSON'S RATIO 

In this section we use the results obtained by Ting and Barnett [4] and relate them to 
regions of the stability triangle. The relation (16) can be written in the form used by Ting 
and Barnett [4] 

_ u( m ,n;s 2 ,s 3 ) _ 2p + Q{s ^ 

e(n;s 2 ,s 3 ) 

withg=[2 (S' u + S' 12 )-S' u }. 

I. Consider a pair of unit vectors: m = (—110) and n = (110). For it, P(m, n) attains 
the maximal value 1/2, and v = Q/A. 
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i) If S' 12 (s2, S3) > 0, Q (s 2 , S3) > 0, then u/e < — (3/4. For s 3 > 0, the inequalities imposed 
on S' 12 and Q define the region deb, whereas for s 3 < - the region aed of ST (cf. Fig. 7). 
Because for all pairs (m, n) Poisson's ratio is negative, elastic materials belonging to these 
two regions are completely auxetics. 

ii) If (s 2 , s 3 ) < 0, Q (s 2 , s 3 ) < 0, then u/e > — Q/A. For s 3 > 0, the inequalities imposed 
on and Q define the region cfe , whereas for S3 < - the region ceh of ST (Fig. 7). 
Because for all pairs (m, n) Poisson's ratio is positive, elastic materials belonging to these 
two regions are non auxetics. 

iii) If S' 12 (s 2 ,s 3 ) > 0, Q(s 2 ,s 3 ) < 0, then u/e < — (5/4. The inequalities imposed on S[ 2 
and Q define the region efb (Fig. 7). Because there exist pairs (m, n) for which Poisson's 
ratio is positive and pairs for which u is negative, elastic materials belonging to these two 
regions are auxetics. 

iv) If S^ 2 (s2,s 3 ) < 0, Q(s 2 ,s 3 ) > 0, then u/e < — Q/A. The inequalities imposed on S' 12 
and Q define the region ache (Fig. 7). Because there exist pairs (m, n) for which Poisson's 
ratio is positive and pairs for which u is negative, elastic materials belonging to these two 
regions are auxetics. 

II. Consider the pair of unit vectors m = (011), n = (100). In this case P(m, n) reaches 
its minimal value 0, and u = — S12, e = sn. 

If S[ 2 (s 2 , s 3 ) > 0, Q (s 2 , s 3 ) > 0, then u < 0. For s 3 < the inequalities imposed on S' 12 and 
Q define the region aed, whereas for s 3 > - the region ebd of ST (Fig. 7). Because for all 
pairs (m, n) Poisson's ratio is negative, elastic materials belonging to these two regions are 
completely auxetics. 

ii) If S' 12 (s 2 , S3) < 0, Q (s 2 , s 3 ) < 0, then u > — si 2 . For s 3 > the inequalities imposed on 
S' 12 and Q define the region cfe, and for s 3 < - ehc (Fig. 7). Because for all pairs (m, n) 
Poisson's ratio is positive, elastic materials belonging to these two regions are non auxetics. 

iii) If S[ 2 (s 2 , S3) > 0, Q (s 2 , S3) < 0, then u/e > — s i2 . For s 3 > the inequalities imposed 
on S[ 2 and Q define the region efb (Fig. 7). Because there exist pairs (m, n) for which 
Poisson's ratio is positive and other pairs for which u is negative, elastic materials belonging 
to these two regions are auxetics. 

iv) If 5(2 (s 2 , s 3 ) < 0, Q (s 2 , s 3 ) > 0, then u/e > s 12 . The inequalities imposed on S[ 2 and Q 
define the region ache (Fig. 7). Because there exist pairs (m, n) for which Poisson's ratio 



11 



is positive and other pairs for which v is negative, elastic materials belonging to these two 
regions are auxetics. 

We conclude that criteria of Ting and Barnett in the unique way divide ST into regions 
having definite auxetic properties. These regions are displayed in Fig. 7. Note that v can 
be isotropic only for complete auxetics. For auxetics Poisson's ratio have to be anisotropic. 



regions-ST.png 



FIG. 7: Auxetic properties of cubic elastic media - characteristic regions of the stability triangle. 
Insert: the enlarged part of ST with some intermediate-valence compounds indicated. Compounds 
distinguished by bold letters are the completely auxetics (cf. [20]). 

Finally, we present results of our calculation of Poisson's ratio, which confirm results 
described in previous sections of our paper. For n =< 110 >, the dependence of v on 
directions m for points belonging to the characteristic regions of ST is presented in Fig. 
8. The solid line represents complete auxetics (the region ebd in Fig. 7). The dotted 
line represents complete nonauxetics (the region c/e), and the broken line - (non-complete) 
auxetics (the region ache). For them Poisson's ratio has to be anisotropic 

Anisotropy properties of Poisson's ratio for several vectors n for the point s = 
(s 2 = —0.2, s 3 = —0.65) G ache region of ST, are shown in Fig. 9. 
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nu-vs-n.png 



FIG. 8: Dependence of u(m, n; S2, S3) on directions of m perpendicular to n = [110] for several 
points belonging to characteristic regions of ST 



nu-on-ST . png 



FIG. 9: Dependence of v on directions of m perpendicular to several directions n for ST point 
S2 = —0.2, S3 = —0.65 belonging to ache - the auxetic region of ST 
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APPENDIX A 

In the Appendix we write down the explicit formulas for Young's, shear modules and 
for Poisson's ratio of cubic media in terms of components of the stiffness and compliance 
tensors. 

1 1 l~3p(n) 1 -p(n) 

E(n) 3(C n + 2C 12 ) 3(C n - C 12 ) + 2C 44 ' 1 ' 

1 P(m,n) | l-2P(m,n) > 



4G(m,n) (C n - C 12 ) 4C 44 

z/(m n) = £(n) ^ + ^t^ 11 ~ Cl2 ~ 2C u)P(™, n ) + ggjjglg (A3) 

2(Cn + 2Ci 2 )(Cn — Ci 2 )C 44 

= 5 12 + i,S 44 + (s n - 5 12 - ^ 44 ) p(n), (A4) 



E(n) 



1 ]s M +(s 11 -S 12 -ls u )P(m,n), (A5) 



4G(m,n) 4 V " 1Z 2 



z/(m, n) = —E(n) 



Si2 + (Sn-S 12 -^Su)P(m,n) 



(A6) 



When S44 = 2 (Cn — C 12 ) or S44 = (Cn — C i2 ) /2 one deals with an isotropic medium. 
If these relations are fulfilled from Eqs. (A1-A6) one obtains expressions for mechanical 
characteristics of isotropic media 

E -i Cn + C i2 g 

(Cn + 2Ci 2 ) (Cn — Ci 2 ) 

C = - (Cn — C12) = 



2 V 11 1Z/ 2(S U -S 12 y 

v = _ES 12 = ° 12 = 

Cn + Ci 2 5*11 
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The knowledge of the components of the compliance tensor of cubic media allows one to 
find the corresponding values of parameters si (I = 1, 2, 3) 

S441S12 + S44S11 — 4S^ 2 + 25'i2<S'ii + 2Sf± 
(25*12 + 5'ii)(5'i2 — S'ii)5'44 

S44S12 + S44S11 + 25*^2 ~~ S12S11 — Sfi 
S44S12 + S44S11 — 45*12 4~ 2<S'i2<S'ii + 2S'ii 

(2S 12 + S n )(2S 12 + ff 44 - 2S n ) 
S44S12 + S44S11 — 45*12 + 2Si 2 Sn + 25ii 
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